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Abstract. We extend and review our analysis of the nucleon —> A(1232) transition electroweak form factors from Ref. [1], 
New fit of the A(1232) vector form factors to electron-proton scattering Fi structure function is introduced as well, leading to 
results different from the popular parametrization of Ref. [2], A clear model dependence of the extracted parameters emerges. 
Fit to neutrino scattering data is performed in all available isospin channels. The resulting axial mass is M aa = 0.85io:os (GeV) 
and C4 (0) = 1 . 10-2: 14. The latter value is in accordance with Goldberger-Treiman relation as long as the deuteron effects are 
included. 
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INTRODUCTION 

The problem of single pion neutrinoproduction (SPP) has been studied for many decades. Its importance has become 
clear with the development of neutrino accelerator experiments, such as MINOS [3], T2K [4], NOvA [5], MiniBooNE 
[ 6 ], and LBNE [7]. In the few-GeV energy region chartacteristic for the above mentioned experiments this interaction 
channel contributes a large fraction of the total cross section. One estimates, that for an isoscalar target and neutrino 
energy of around 1 GeV SPP accounts for about 1/3 of the interactions. 

These SPP events give rise to the background in measurements of quasi-elastic neutrino scattering off nuclear 
targets if subsequent pion absorption occurs. For experiments aiming at electron neutrino appearance measurement 
neutral current K° production process adds to the background in water Cherenkov detectors. Correct understanding 
and modeling of the cross-sections for the SPP is crucial for precise extraction of neutrino oscillation parameters in 
long baseline experiments. 

Theoretical modelling of the SPP processes on nuclear targets is biased by systematic errors coming from nuclear 
model uncertainties. They are driven by the strong nature of hadron interactions inside the nucleus, which do not 
allow for a feasable, exact solution of the problem. Experimental measurements also suffer from these efects. An 
apparent tension between the MiniBooNE and very recent MINERvA SPP data on (mostly) carbon target (Ref. [ 8 , 9]) 
is one of the key examples. For the purpose of analysis of Nucleon (N) to A(1232) resonance transition vertex one 
desires measurements of the neutrino-production on free or almost free targets. Such data exist only for ~30 years old 
Argonne National Laboratory (ANL) [10, 11] and Brookhaven National Laboratory (BNL) [12, 13] bubble chamber 
experiments, where deuteron and hydrogen targets were utilized. In this case one may hope to reduce the many-body 
bias in a reasonable manner with a simple theoretical ansatz [14]. 

One can not understand the neutrino SPP data without introducing an appropriate nonresonant background, see 
Ref. [15]. More recent studies of weak SPP fit the N —> A transition axial form factors utilizing only the neutrino- 
proton channel + p —> p~ ± n + + p [16, 17, 18, 19, 20, 21], Simple total cross section ratio analysis shows, that 
the background contribution is much larger in neutrino-neutron channels. The neutrino-proton SPP channel can be 
described well within a model that contains the A(1232) resonance contribution only, see e. g. Ref. [22]. Thus a 
quantitative, statistical, validation of any pion neutrinoproduction model should be done using all available isospin 
channels, where one gets more information about the background contribution. In Ref. [22] a consistent fit for both 
ANL and BNL data sets with deuteron effects included yielded Cj(0) = 1.19 ±0.08 and M A = 0.94 ±0.03 GeV. 
The attempt to extract the leading Cf (Q 2 ) N —t A form factor parameters in a model containing both nonresonant 
background and deuteron effects has been done in Ref. [19]. The results gave the value of C 5 (0) = 1.00 ±0.11 far 
from the Goldberger-Treiman relation estimate of C 5 (0) ~ 1.15 [23, 24], Lack of the deuteron effects pushes the fit 
even further away from theoretical prediction, giving Cj (0) = 0.867 ± 0.075 in Ref. [17]. From the above mentioned 


models only those in Refs. [16, 18] have been directly tested against the electroproduction data. Authors of [17, 19, 20] 
use vector form factor parametrization from Ref. [2], based on the MAID analysis [25]. The A form factors have been 
fitted directly to the A helicity amplitudes from MAID. Thus the approach from Ref. [2] relied only on the A resonance 
excitation without any nonresonant background. The problem is that the resonance helicity amplitudes extraction 
procedure is model-dependent. One has to make some assumptions how to separate A and background contributions 
from the data and A - background interference effects are strong. We show later that if one uses form-factors obtained 
within one description of pion electroproduction in a model with different physical components the resulting cross 
sections may become imprecise. 

The above mentioned caveats of previous analyses have motivated us to propose an improved approach. We adapt 
and develop the statistical framework of Ref. [22] in order to fit both vector and axial form factors of the A(1232) 
resonance. We use inclusive electron-proton scattering data for the electromagnetic interaction in the A(1232) region, 
adding a new fit of the form factors from Ref. [2]. The deuteron bubble chamber data of ANL and BNL experiments 
are used for the weak interaction. For this analysis we are the first authors to incorporate the neutron channels. In this 
manner we include the data sets, that are very sensitive to the nonresonant background. 


FORMALISM 

In the neutrino single pion production off free nucleon targets one distinguishes three isospin channels: 

v n(l)+ p(p) ~M (0 + n + (k) +p(p') ( 1 ) 

V/j(/)+h(p) —> p (l r )+ n°(k) + p(p') ( 2 ) 

v p(l)+ n ( P ) —> p ( l') + n + (k) + n(p' ) (3) 

with /, p, p' and k being the neutrino, muon, initial nucleon, final nucleon and pion four momenta respectively. The 
definition of four momentum transfer is following: 

q = l — l' = p' + k — p, Q 2 = —q 2 , = (<?°,q) (4) 

and the square of hadronic invariant mass is: 

W 2 = (p + q) 2 =(p' + k) 2 . (5) 

Metric convention g^ v = diag(+, —, —, —) is used throughout this paper. 

For the pion electroproduction we are interested in proton target reactions; 

e~(l )+p{p) —> e~(l') + K + {k) +n(p') (6) 

e~{l)+p(p) -> e~(l') + n°(k)+p{p'). (7) 

The proton interaction channel (1) is dominated by the intermediate A ++ resonance excitation, which makes it very 
sensitive to the properties of this resonance. Neutron channels (Eqs. (2) and (3)) contain a large contribution of 
nonresonant pion production, thus they present more challenges for theorists. They are crucial to verification of any 
consistent SPP model. 


N — y A(1232) transition 

We treat the A(1232) resonance excitation within the isobar framework. The most general form of positive parity 


spin -1 particles electroweak excitation vertex can be expressed as: 
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A relevant information about the inner structure of the A(1232) resonance is contained in a set of vector and axial 
form factors cj’. In this paper they are assumed to be functions of Q 2 only (with the exception of C'J which depends 
also on W). 


Vector contribution 


The conserved vector current (CVC) hypothesis gives the relation between weak and electromagnetic vector form 
factors. In the case of A(1232) resonance and hereby used convention both sets are exactly the same. The size and 
excellent accuracy of the electromagnetic data set allows for an introduction of multiple fit parameters. 

We explore two models of vector form factors. The first parametrization, refered to as “Model I”, has the same 
functional form, as in Ref. [2]: 
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In the above equations My = 0.84 GeV is the standard vector mass. Everything else is treated as a fit parameter. 

We also propose our own model of electromagnetic form factors. We assume that the N —> A transition form 
factors have the same large Q 2 behaviour as the electromagnetic elastic nucleon form factors. There exist theoretical 
arguments [26] suggesting that at Q 2 —> °° the nucleon form factors fall down as 1 / Q 4 . Following these assumptions 
we adopt appropriate Pade type parametrization used previously to parametrize the electromagnetic form factors 
of the nucleon [27]. In this manner we allow for a deviation from the St/(6)-symmetry quark model relations 
C\ (Q 2 ) = —( M/W)C 3 ( Q 2 ) and Cj = 0 between the form factors [28]. Finally, we assume the dipole representation 
of C' (Q 2 ) to reduce the number of parameters. Altogether, our parametrization has the following form: 
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This parametrization reproduces quark model relation between C) and C\ at Q 2 = 0. It is also allows for nonzero 
value of S \/2 helicity amplitude. We call it “Model II”. 


Axial contribution 

Here the leading contribution comes from C~ (Q 2 ) which is an analogue of the isovector nucleon axial form factor. 
Partially conserved axial current (PCAC) hypothesis relates the value of C) (0) with the strong coupling constant f* 
through off-diagonal Goldberger-Treiman relation [23, 24]: 
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but we will treat Cj (0) as a free parameter. We assume, that Cf has a dipole Q 2 dependence: 


C5 (G 2 ) 


cf( 0 ) 



The axial mass parameter Maa can be related to “resonance axial charge radius”. It is also subject to fit, but we expect 
it to be of the order of 1 GeV. 

The Cg form factor is an analogue of the nucleon induced pseudoscalar form factor. One can use PCAC to relate it 

toC^: 


c e(Q 2 ) 


M 2 

ml + Q 2 


Cf(Q 2 ), 


(19) 


where m n is average pion mass. The C^(Q 2 ) is the axial counterpart of the very small electric quadrupole (E2) 
transition form factor. Unfortunately, bubble chamber data set is too inacurate to precisely measure its effect. Due 
to expected similarities between A and nucleon properties we set C 3 = 0. For the C 4 we use the Adler model relation 
[29]: 


C^(G 2 ) = -C^(G 2 )/4. (20) 

In this way the axial contribution is fully determined by C^(Q 2 ). Altogether there are two free parameters: C 5 (0) and 
Maa- If there were enough experimental data one could drop the Adler relation and treat C 4 (Q 2 ) as an independent 
form factor or even determine, whether Cf(Q 2 ) has nonzero value. However, the ANL and BNL experimental data do 
not have sufficient statistics even to obtain separate fits of Cf and C 4 [30], see also the discussion in Ref. [19]. 
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FIGURE 1. (Color online) Basic pion production diagrams from [17]: a) Delta pole (AP), b) crossed Delta pole (CAP), c) nucleon 
pole (NP), d) crossed nucleon pole (CNP), e) pion-in-flight (PIF), f) contact term (CT) and g) pion pole (PP). 

We express the inclusive double differential SPP cross section for neutrino scattering off nucleons at rest as: 


d~a 


1 


jGj. cos 2 Gc 


W r 

M 2 J ' 


d 3 k 


dQ 2 dW 32 kE 2 f '~‘ j M 2 J (2n) 3 2E K {k)E(p’) 

L pv = / P / /v +/ v zV _^v / . z / + - e pva/3^^ 

A^ v = Y.(K N, \j?M\N)(7tN'\j v cc mN)*- 


jrLpvA FV 8(E(p') +E K (k) —M — q°) 


spins 


( 21 ) 


























where E is the incident neutrino energy, M is the averaged nucleon mass, E n {k) and E(p') are the final state pion and 
nucleon energies, Gf = 1.1664 -10 11 IVIeV 2 is the Fermi constant, G' v - the leptonic and A“ v - the hadronic tensors. 
The Cabibbo angle, cos (6c) = 0.974, was factored out of the weak charged current definition. 

The information about dynamics of SPP is contained in transition matrix elements, (nN' |y'£,(0)| N), between an 
initial nucleon state \N) and a final nucleon-pion state | kN'). 

In the model of this paper the dynamics of SPP process is defined by a set of Feynman diagrams (Fig. 1) with 
vertices determined by the effective chiral field theory. They are discussed in Ref. [17], where one can find exact 
expressions for j^ 1 . The same set of diagrams describes also pion electroproduction, with the exception of the pion 
pole diagram, which is purely axial. We call this approach "HNV model" after the names of the authors of Ref. [17]. 


Deuteron effects 


In this paper we consider a deuteron model based on phenomenological nucleon momentum distribution, /(/?), 
taken from the Paris potential [31]. It is assumed that the spectator nucleon does not participate in the interaction 
and that there are no final state interactions (FSI). This assumption is based on the results of Ref. [32] where for the 
quasielastic scattering case FSI were shown to be negiligible as long as the neutrino energy is larger than 500 MeV. 
The very recent study (newer, than the hereby analysis) of Ref. [33] proved the FSI to be important in the mt + channel, 
where they lead to a substantial reduction of the cross section for forward-going pions. Further studies regarding the 
impact of FSI effects are needed, also in the ANL and BNL experimental data analysis where the event selection is 
based on spectator approach as well. We introduce also the effective binding energy: 

B{p) = 2E{p)-M D , (22) 


where Mu is deuteron mass. The expression for the cross section becomes: 
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with qf 1 = ( q° —B(p),q ) and v re i = \J (/ • p) 2 /EE(p). We also define a Jacobian \ 
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whose explicit form is complicated because the invariant mass W depends both on the energy transfer q° and the lepton 
scattering angle ©. 


RESULTS OF THE ANALYSIS 

Our main goal is to have a reliable model of weak pion production. Because the neutrino SPP data are sufficient only 
to obtain information about leading axial coupling of the A(1232) resonance, we assume that the extraction of vector 
and axial form factors can be done independently using first respective electron scattering and then neutrino SPP data. 
In the next paragraphs we describe details of our procedure. 


Vector Contribution to Weak SPP and Electroproduction 

Our aim is (due to a poor quality of the neutrino SPP data) to reproduce correctly only the most important 
characteristics of the neutrino SPP reactions. These include overall cross sections and distributions in Q 2 . Detailed 
analysis of the electroproduction data bases on pion angular distributions. Such a task is beyond the scope of this 
paper. 

We use the information contained in electron-proton Fi data from [34]. We include 37 separate series of /•? data 
points from the lowest value of Q 2 (0.225 GeV 2 ) up to 2.025 GeV 2 . This Q 2 range overlaps with the one in ANL data. 



The data desribe the inclusive structure function, thus we limit our fit to values of invariant mass W up to M p + 2m n . 
Beyond that value the experimental data contain two pion production and then more inelastic channels. With this 
limitation for Q 2 < 2.025 and W < M p + 2m K we are still left with 603 data points. 

In order to ensure that the results will reproduce well the data at the A( 1232) peak we expanded our fit to higher value 
of invariant mass W = 1.27 GeV. There are no exclusive electron SPP data in the region W £ (M + 2m Jl . 1.27 GeV). 
Thus we chose to add to our fit a term in which MAID 2007 model predictions are taken as 228 fake data points 
with errors identical to respective Osipenko et al. [34] points. We could not apply the MAID model directly in our 
fits since the exact formulas for their SPP amplitudes have never been published. These additional points have been 
generated using the on-line version of MAID (http : / /wwwkph . kph . uni-mainz . de/MAID/ /). We have also 
included the information about MAID 2007 model helicity amplitudes. The caveat is that the pion electroproduction 
experimental results contain both resonant and nonresonant contributions (see e.g. Ref. [35] ). Thus the extracted 
helicity amplitudes depend on how one defines the "Delta" and "background". The HNV model differs with MAID 
in the treatment of both. One can not expect the resulting helicity amplitudes to be the same. From that reason the 
information about helicity amplitudes has been given a large ad hoc error assumption in our estimator. 


Results 

The best fit results of our vector form factor parametrization given by Eqs. (11-13) and Eqs. (14-16) are shown 
in Table 1. We also present there the values from Ref. [2] in order to compare directly with our model I. In both 
models the best fit value of C) (0) is close to the one from Ref. [2] and we get a clear beyond-dipole Q 2 dependence 
of Cj(0 2 ) and C'J (Q 2 ). Surprisingly, the Q 2 dependence of C y (Q 2 ) is exactly dipole (1 +Q 2 /My )~ 2 in model II. 
Most importantly, we have shown, that extracted form factors are model-dependent. This clearly follows from the 
difference of best-fit parameters between our model I and their counterpart from Ref. [2]. One can see that, besides 
the similarity in the leading form factor value C' y (0), both fits differ by large. 


TABLE 1. Best fit coefficients for vector form factors given by Eqs. (11-13) (“Model I”) 
and Eqs. (14-16) (“Model IF’). We do not report let errors because of hybrid character of 
our estimator, see explanations in the text. 
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Ref. [2] 
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MODEL I 

2.00 

-6.77 

0.68 
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MODEL II 
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4.73 

-0.39 
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0.13 

1.68 
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Fig. 2 shows that qualitatively in the region below two pion production threshold our fit reproduces the data rather 
well. Our form factors lead to better agreement with the electron scattering data than the form factors considered 
in Ref. [2], The same trend is clearly seen in Fig. 3, where our best fit results are compared to the inclusive electron- 
proton scattering cross section data. Inspection of Fig. 2 shows that biggest disagreement with data is exhibited in 
region of low W. Our fits are going to be used in the analysis of neutrino scattering data and some discrepancy at low 
W is of no practical importance. 

Another conclusion is that each physical model of single pion production needs its own separate resonance form 
factor analysis. Any change of description of one of the elements such as A propagator and width, background 
amplitudes, unitarity constraint etc. will affect the results. In other words, the HNV model should be used together 
with vector form factors fited using the HNV model in order to increase the accuracy of its predictions. 

Because both proposed form factor sets lead to very similar results, we choose to use “Model II” in the axial fits. 


Axial Contribution to weak SPP: fits to bubble chamber data 

We consider a statistical framework, proposed in Ref. [22], ANL used a neutrino beam with mean energy below 
1 GeV and a large flux normalization uncertainty A panl ~ 20% that was not included in the published do/dQ 2 cross 
section for the reaction in Eq. (1) [22]. ANL reported the data with the invariant mass cut W < 1.4 GeV, which allows 
us to confine to the A(1232) region. We can neglect contributions from heavier resonances, whose axial couplings are 
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FIGURE 2. (Color online) Best fit results for vector form factors given by Eqs. (11-13) (“Model I”) and Eqs. (14-16) (“Model 
II”) plotted against experimental data from Ref. [34] and predictions of HNV model with original Lalakulich-Paschos form factors 
of Ref. [2], Vertical lines show the 2 n production threshold. 


by large unknown. Our analysis uses information from all available SPP isospin channels. The detiled description of 
the statistical approach can be found in Ref. [1]. 

We treat Cj (0), M AA and normalization factor panl as free fit parameters in the analysis of axial N —> A(1232) 
transition. We present our results in Tab. 2. where fits to all three channels separately as well as the joint fit to three 


TABLE 2. Best fit for the A( 1232) axial form factors on deuteron target. Errors 
for C 5 (0) and M A a were obtained after marginalization of Panl- 


Channel 

C£( 0) M AA ( GeV) p ANL X 2 /NDF NDF 

Vp+P~> +P + n + 

Vp +n -»• p~+p + K° 
Vfi -{- n —^ p -f- n -\- 7l^ 

Joint fit 

1.1 l±o:15 0.97±iu7 1.04 0.20 6 

1.31 ±0:77 1.00±8:S 0.93 1.52 9 

2.83±a«s 0-76T8:!! 0.94 1.47 9 

1.10^14 0.85±8:81 0.90 2.06 30 


channels are listed. In each case the number of degree of freedom is calculated as: 

NDF = No. Q 2 bins — No. fitted parameters. 

We see, that taken separately the pn + (Al) and pn° (A2) channels are statistically consistent, albeit their predicted 
scale parameters differ by around 10%. The latter channel seems to carry less information on the N —> A transition 


























































FIGURE 3. (Color online) Comparison of our best fit results and HNV model with Lalakulich-Paschos form factors of Ref. [2] 
plotted against inclusive p(e,e') data (not included in the fit) from Ref. [36] (left panel) and Ref. [37] (right panel). The Q 2 values 
at peak are from left to right: 0.1 (GeV 2 ) and 0.95 (GeV 2 ) respectively. 


axial current than the first one. This fact is reflected in larger uncertainties. We explain it by a bigger background 
contribution to that channel, which makes it less sensitive to changes in the A resonance form factors. 

The mi + (A3) channel gives results inconsistent with the other two. Cf (0) is obtained twice as large as for the pn 
and p7t° channels and Maa significantly smaller. Here the number of events reported by ANL is comparable to pK {) 
channel, but theoretical cross section predicted by our model are smaller. This results in the overestimation of Cj (0). 
Surprisingly, the fits to separate isospin channels give acceptable values of Xmin f° r both neutron channels. 

Deuteron effects affect mostly the value of Cj (0) by up to 20%. In the joint fit on free proton and neutron targets we 
obtained Cf( 0) = 0.93 4 S:n and Maa = 0.8 1 '<!:!« GeV compared to 1.10 and 0.85 GeV in the deuteron case. A significant 
improvement with respect to previous fits to HNV model done in Refs. [17, 19] is that with deuteron target effects we 
get the best fit value of Cf( 0) within 1 a range from the theoretical Goldberger-Treiman relation. The joint fit agrees 
also on the lcr level with separate fits on pn {) and pn + channels. 



FIGURE 4. (Color online) Differential cross sections on a free nucleon target for our best fit form factors and for fits of Ref. [19]. 

We have compared our best fit form factors with the results of Ref. [19] for two different neutrino energies 
characteristic for Minos and MINERvA experiments in Figs. 4 and Fig. 5. Our form factors predict different shape 
and magnitude of do/dQ 2 . The size of effect is bigger for the lower values uf Q 2 and for proton channel, where the A 
contribution dominates. In the latter one has up to 10% difference in the cross section. 

Normalization factors panl fitted separately for each channel are different for neutrons and protons. The proton 
channel prefers the data to be scaled up and both neutron channels prefer the data to be scaled down. The joint fit uses 
the same panl parameter for all channels and seems to prefer the data to be scaled down even more ( panl ~ 0.90). 
The values of panl are all within the assumed flux normalizadion error A panl- 
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FIGURE 5. (Color online) Same as Fig. 4 but for E v = 3 GeV. 


Finally, we noticed that the best fit values for C 5 (0) and Maa in the pn + (Al) channel are different from those 
obtained in Ref. [22]. We explain it by the inclusion of nonrenonant background in the current analysis. 


CONCLUSIONS 

In this paper we review a new attempt to get an information about weak N —> A transition matrix elements previously 
presented in Ref. [1] as well as we extent this analysis by presenting the newest fits of the vector form factors and 
discussion of the axial form factors. The fit to electromagnetic has clearly shown, that the extracted vector form 
factors of the A resonance are model-dependent, e. g. the HNV model gives the best results with the form factors 
extracted using the full HNV model, as it has been done in this paper. 

We discussed axial form factor fits obtained based on the analysis of all three neutrino SPP channels, including the 
neutrino-neutron channels. In previous works usually only neutrino-proton channel was utilized to extract the axial 
form factors. A critical analysis of neutrino-neutron channel, on qualitative level, appears also in other papers see 
e. g. [15, 17], a detailed discussion can be found in Ref. [20] as well. The obtained value of Cj (0) agrees, on the la 
level, with the Goldberger-Treiman relation if the deuteron effects are taken into account, which is an important result. 
If one neglects the nuclear effects the resulting C 5 (0) value is lower. Also, there is a strong tension between nil ' 
and remaining two channels (see also [17] and [20]). The same theoretical model does not seem to give a consistent 
reproduction of data in all channels. 

There can be various reasons for that. Firstly the existing bubble chamber data on neutron SPP channels are of poor 
statistics. Secondly, the chiral model for the background is well justified only near the pion production threshold. It 
may be not reliable in the A(1232) peak region. Last, but not least, the iik 1 channel is subject to large FSI effects, 
as shown in Ref. [38]. Thus the spectator model used both in experimental analyses of ANL and BNL as well as in 
our calculations may give invalid results in this channel. The plots in Ref. [38] suggest a reduction of the n7t + cross 
section due to FSI. Further studies are needed. 

Still another reason of these difficulties may come from a missing unitarization of the model. This constraint, 
following the Watson theorem [39], imposes a relation between phases in lepton-nucleon and pion-nucleon elastic 
scattering amplitudes. Unitarity is not satisfied in our approach. In a recent study Nieves, Alvarez-Ruso, Hernandez 
and Vicente-Vacas [40] tried to correct the HNV model by introducing phenomenological phases. They obtained 
better agreement of the best fit value of Cf (0) with the Goldberger-Treiman relation 1 . This is a strong indication of 
the importance of proper model unitarization for the pion neutrinoproductioncase. More studies are necessary. 

Relatively large values of experimental errors and our inability to extract independent C 3 and C 4 form factors 
implies, that better statistics SPP measurements in the A region on proton or deuteron targets are badly needed. 
Keeping in mind difficulties in the treatment of nuclear effects on heavier targets it is the only way to get precise 
information about the N —> A axial transition matrix elements. 


See also the contribution of E. Hernandez presented during the CETUP* 2014 workshop. 
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